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ABSTRACT 



Aims. The Galileon model is a modified gravity theory that may provide an explanation for the accelerated expansion 
of the Universe. This model does not suffer from instabilities or ghost problems (normally associated with higher-order 
derivative theories), restores local General Relativity - thanks to the Vainshtein screening effect - and predicts late 
time acceleration of the expansion. 

Methods. In this paper, we derive a new definition of the Galileon parameters that allows us to avoid having to choose 
initial conditions for the Galileon field, and then test this model against precise measurements of the cosmological 
distances and the rate of growth of cosmic structures. 

Results. We observe a small tension between the constraints set by growth data and those from distances. However, 
we find that the Galileon model remains consistent with current observations and is still competitive with the ACDM 
model, contrary to what was concluded in recent publications. 

Key words. Supernovae: general - Cosmology: observations - Cosmology: dark energy 
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1. Introduction 

The discovery of th e accelerated expansion o f the Universe 
(|Riess et al.l . 119981 : iPerlmutter et all I1999D led cosmolo- 
gists to introduce dark energy to explain our Universe. 
Adding a cosmological constant (A) to Einstein's General 
Relativity is the simplest way to interpret observational 
data. However, even if adding a new fundamental constant 
is satisfactory, the value of A obtained from numerous mea- 
surements results in significant fine tuning and coincidence 
issues. Thus, there is theoretical motivation to find alter- 
native explanations, such as modified gravity models. 

The Galileon model is just such a formulation. It was 
first proposed bv lNicolis et al. ( 2009) as a general theory in- 
volving a scalar field, hereafter called ir, and a second order 
equation of motion invariant under a Galilean shift symme- 
try (d^n — > dfj_ir + 6 M , where is a constant vector). This 
symmetry was firs t noticed in braneworld theories such as 
the DGP model of lDvali et all (|2000h . The DGP model has 
the advantage of providing a self-accelerating solution to 
explain the expansion of the Universe, but which is plagued 
by ghost and instability issues. Galileon theories are a gen- 
eralization of the DGP model which avoids these problems. 
The G alileon model was then derived in a covariant formal- 
ism bv iDeffavet et alJ (|2009h . It was also shown that this 
model forms a subclass of the general tensor-scalar theo- 
ries involving only u p to second order derivatives originally 
found by Horndeski (Horndcski. I1974T ). 



In a 4-dimension spacetime, only five Lagrangian terms 
are possible when forming an equation of motion for ir 
invariant under the Galilean symmetry. Therefore, the 
Galileon Lagrangian has only five parameters. In the 
Galileon theory, as in the DGP th eory, a scre e ning mecha- 
nism called the Vainshtein effect (|Vainshteinl . Il972f ) arises 
near massive objects due to non-linear derivative self- 
couplings of the 7r field. These ensure that the Galileon 
fifth-force is screened near massive objects, and preserves 
General Relativity on local scales where it has been exper- 
imentally tested to high precision. However, this screening 
is only effective below a certain distance from massive ob- 
jects (the Vainshtein radius) which depends on the mass 
of the object and on the values of the Galileon parameters 
([Burrage fc Seervl , |201Q[ ). Experimental constraints on the 
Galileon paramete rs based on local tests of gravity have 
been proposed by IBrax et all (|201ll ) and IBabichev et al.l 

(poll . 

Recently, the Galileon model has been tested aga i nst ob - 
servational c o smolo gical data by Appleby fc Linderl ((2012), 
lOkada et all (|2012j ). and lNesseris et alJ (|2010l) . These au- 
thors tend to reject the Galileon model due to tensions be- 
tween growth of structure constraints and the other cosmo- 
logical probes. The evolution of the Universe in the Galileon 
theory is based on differential equations involving the n 
field, which requires one to set initial conditions, and the 
above studies resorted to different methods for setting these 
initial conditions. In this work, we avoid this issue by intro- 
ducing a new parametrisation of the Galileon model which 
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renders it independent of initial c onditions. Combined with 
the oretical constraints derived in [Appleby fc Lindeil (|2012t ) 
and iDe Felice fc Tsuiikawal f|2011h , we compare our model 
with cosmological observables, and find that the Galileon 
model is not significantly disfavoured by current observa- 
tions. 

In this paper, we use the most recent measurements 
of SN la luminosity distances, the cosmic microwave back- 
ground (CMB), and baryon acoustic oscillations (BAO). 
The highest quality SN la sample currently a vailable is the 
SNLS 3 sample described inlGuv etaD (|2010ft : IConlev et al.l 
(|201ll) ; ISullivan et all (1201 ll). For the CMB we u se the ob- 
servables from WMAP7 ~(|Komatsu et all . 1201 ll). and the 
set of BAO distances of the BOSS analysis ([Sanchez et al.l . 
l2012t) . The growth of structures is an important probe for 
distinguishing modified gravity models from standard cos- 
mological models such as ACDM, so it has to be used care- 
fully. In this work, we use fas(z) measurements from sev- 
eral surveys, corrected from the Alcock-Paczynski effect. 

Section [5] provides the Galileon equations used to com- 
pute the evolution of the Universe and the theoretical con- 
straints imposed to the Galileon field. Section [3] describes 
the likelihood analysis, data samples and how cosmological 
observables are computed. Section U gives the constraints 
on the Galileon model derived from data, and Sect. [S] dis- 
cusses these results and their implications. We conclude in 
Sect. El 

I 



2. Cosmology with Galileons 

2.1. Lagrangians 

The Galileon model is constructed from the assumption 
that the scalar field equation of motion is invariant un- 
der Galilean symmetries: d^ir — > d^ir + 6 P , where b p 
is a constan t four vector. By imposing this symmetry, 
iNicolis et al.l (|2009f) showed that there are only five pos- 
sible Lagrangian terms Li for the Galileon model ac- 
tion. The covari ant formulation of the Galileon Lagrangian 
was derived in iDeffavet et al.l (|2009D . In this paper we 
st art with this covar i ant a ction with the parametrisation 
of lApplebv fc Linderl IpOll : 



S 




1 5 

IE 



CiLi 



(1) 



with L m the standard matter Lagrangian, Mp the Planck 
Mass, R the Ricci scalar and g the determinant of the met- 
ric. The CjS are the arbitrary dimensionless parameters of 
the Galileon model, which weight the different terms. The 
G alileon Lagrang i ans h ave a covariant formulation derived 
in IDeffavet et all (|2009t ): 



L x = A/ 3 7r, L 2 = (V M 7r)(V"7r), L 3 = (□ 7 r)(V A1 7r)(V^^)/Af 3 
La = (V, l7 r)(V^7r) [2(Oir) 2 - 2ti,^ v - fl(V M 7r)(V*7r)/2] /M 6 
L 5 = (V^XV'Tr) [(Dtt) 3 - 3(Dtt)tt,^tt^ + 2^ t >»ir ;v ^tt, p ^ - 6tt,^^ p G up ] /M 9 



(2) 



where M is a mass parameter defined as M 3 = H^Mp, 
where Hq is the current value of the Hubble parameter. 
With this definition the CjS are dimensionless. 

Z/2 is the usual kinetic term for a scalar field, whereas 
L 3 to L$ are non-linear couplings of the Galileon field to 
itself, to the Ricci scalar R and the Einstein tensor 
providing the necessary features to modify gravity and 
mimic dark energy. L\ is a tadpole term which acts as the 
usual cosmological constant, and may furthermore lead to 
vacuum instability as it is an unbound ed potential term. 
There fore, in the following we set c\ = 0. Appleby fc Linderl 
(2012) propose additional direct linear couplings to mat- 
ter to add to the action: a linear coupling to matter 
Lq = cqttT^ /Mp and a derivative coupling to matter 
Lq = coduirdyirT^ / ( MpM 3 ) which arises in some bran e- 
world theories (see e.g. iTrodden k, Hinterbichlerl ([201 lM . 
In this work, we do not consider these optional extensions 
to the theory, as they would add more free parameters cq 
and cg to the Galileon model. 

This action leads to three differential equations : two 
Einstein equations ((00) temporal component and (ij) spa- 
tial component) coming from the variation of the action 
with respect to the metric g pv , and the scalar field equation 
of motion from the variation of the action with respect to 
the 7T field. The equations are given explicitly in Appendix 
B of lApplebv fc Linderl (|2012l ). With these three differential 



equations the evolution of the Universe and the dynamics 
of the field can be computed. 

To solve the cosmological equations, we choose the 
Fricdmann-Lemaitre-Robertson- Walker (FLRW) metric. 
With no direct couplings, the functions to compute are the 
Hubble parameter H = a/a (with a the cosmic scale fac- 
tor), and x = it' /Mp, w ith a prime denoting d/dlna (see 
lApplebv fc Linderl (|2012h and Sect. HOI). 

2.2. Initial conditions 

To compute the solutions of the above equations, we need 
to set one initial condition for x. We arbitrarily choose 
to define this initial condition at z = 0, which we denote 
xq = x(z = 0). Unfortunately, we have no prior information 
about the value of the Galileon field nor its derivative at 
any epoch. Fortunately, xo can be absorbed by redefining 
the CiS as follows: 

Cj = CiXg (3) 

x = x/xq (4) 
H = H/H . (5) 
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This redefinition allows us to avoid treating xq as an ex- 
tra free parameter of the modeQ- Doing so, the qs remain 
dimensionless, and the initial conditions are simple: 



x 



1, H = 1. 



(6) 



Note that the (00) Einstein equation could also be used 
as a constraint equation to fix xq (see Appendix A) given 
a set of cosmological parameters as, an d Doing 
so, we would observe a degeneracy between the parame- 
ters: the same cosmological evolution can be obtained with 
small as and a high xq, or with high CjS and a small xq. In 
other words, different sets of parameters {ci,xo} produce 
the same cosmology, i.e. the same p n {z), which is undesir- 
able. Our parametrization avoids this issue by absorbing 
the degeneracy between the c^s and xq into our qs. 



with iVe g = 3.04 the standard effective number of neutrino 
species (jManeano et all I2002D . The photon energy density 



at the current epoch is given by f2°/i 2 = 2.469 x 10 
(where, as usual, h = -Ho/(100 km/s/Mpc) for X"cmb = 
2.725 K. 



2.4. Perturbation equations 

To test the Galileon model predictions for the growth 
of structures, we also need the equations describ- 
ing density pertu r bation s. We follow the approach of 
Appl eby fc Lindei] (120121) for the scalar perturbation. 
Appleby fc Linden ((201 2T ) performed their computation in 
the frame of the Newtonian gauge, for scalar modes in the 
subhorizon limit, with the following perturbed metric: 



2.3. Cosmological equations 

To compute cosmological evolution in the Galileon model, 
for simplicity we assume that the Universe is spatially flat, 
in agreement with current observations. In this work, we 
use the FLRW metric in a flat space: 



ds 2 = -dt 2 



a 2 5jjdx l dx-' . 



(7) 



When writing the cosmological equations, we can mix the 
(ij) Einstein equation and the it equation of motion to ob- 
tain the following system of differential equations for x and 
H: 



H' 



aX — (77 
+ 

ap — acu 
W7 — A/3 
a/3 — auj 



with 



a = — Hx - 3c 3 H 3 x 2 + 15c 4 H 5 x 3 - —c 5 H 7 x 4 
6 2 



7 - ^-H 2 x - c 3 H 4 x 2 + ^-c 5 H s x 4 



(-2 



H 2 - 2c 3 H*x + 9c 4 H°x 2 - 10c 5 H a x 



8-3 



(8) 
(9) 

(10) 

(11) 
(12) 



a = 2H + 2c 3 H 3 x 3 - 15c 4 H 5 x 4 + 2lc 5 H 7 x 5 (13) 

(14) 



X = 3H 2 + ^ + %H 2 x 2 - 2c 3 H 4 x 3 



y c 4 ff 6 x 4 - 9c 5 H s x 5 

6^3 



u = 2c 3 H*x 2 - \2c 4 H x 6 + 15c 5 H s x 



(15) 



as derived in the formalism of lApplebv fc Linderl (|2012l) . 
but using our normalisation for the qs. We obtain the 
same equations except that the qs are changed into c^s, 
and that we have a different treatment for the initial con- 
ditions. Equations [8] and [9] depend only on the CjS and 
f2°. The radiation energy density in equation 1141 is com- 
puted from the usual formula = + 0.2271 AW) 



1 If the optional coupling parameters Co and cg are included, 
they should be redefined as: cq = cqXq and cq = cgXq. But if 
Co 0, two initial conditions are needed (no and n' ). With our 
new parametrization we would have to introduce and fit a new 
parameter ro = ho/tt'q. 



ds 2 = -(1 + 2i>)dt 2 + a 2 (l - 2cj))5ijdx l dx j (16) 

In this context, the perturbed equations of the 
(00) Einstein equation, the (ij) Einstein equation, the 7r 
equation of motion and the equation of state of matter, in 
the quasistatic approximation, are, respectively: 



1 



K 4 \7 2 tfj - n 3 \7 2 (j) = Kl \7 2 6y 



K5V 8y — K4V 



1 



H 2 M 2 



-K5V if) — KiV 4> — K 6V Sy 



(17) 
(18) 
(19) 



H 2 C + HH'S' m +2H 2 8' m = (20) 

where Sy — Sn/Mp is the perturbed Galileon, V = V/Hq, 
p m is the matter density, S m = 5p m / p m is the contrast 
matte r density, the same as in lApplebv &; Linderl 

(|2012f ). but rewritten following our parametrization: 



ki = -6c 4 H 3 x 3 H'x + Hx 



Hx 



(21) 



«3 



c 5 H 5 x 3 (12Hx' + IhH'x + 3Hx) 

-l-^H 4 x 4 -3c 5 H 5 x 4 {H'x + Hx') (22) 

(23) 
(24) 



k 4 = -2 + 3c 4 iPx 4 - 6c 5 H°x 



6^5 



K 5 = 2c 3 H 2 x 2 - \2c 4 H*x A + 15c 5 H°x 

n 6 = y - 2c 3 (H 2 x' + HH'x + 2H 2 x) 

+ c 4 (\2H 4 xx' + 18H 3 x 2 H' + 13H 4 x 2 ) ( 25 ) 
- c 5 (18H 6 x 2 x' + 30H 5 x 3 H' + 12H 6 x 3 ). 

With equations [T7] to [201 we can obtain a Poisson equa- 
tion for tp, with an effective gravitational coupling 
which varies with time and depends of the Galileon model 
parameters CiS: 



C 



eff 



■"0 



4(k 3 k 6 - k 2 ) 



K5(k 4 Ki — K5K3) — K 4 (k 4 Kq — K5K1) 



(26) 
G N (27) 
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with Gm Newton's gravitational constant. These equations 
can be used to compute the growth of matter perturba- 
tions in the frame of the Galileon model (see Sect. I3.2.4|) . 
Tensorial perturbations modes also exist, and are studied 
in Sect. I2TQ1 



2. a Laplace stability condition for the propagation speed 
of the perturbed field: 

2 _ 4kik 4 k 5 - 2k 3 kI - 2k 2 a k & 

c s — 77, I o 2\ ^ U ^ 6 I 

K4(2K4K2 + 0K5J 



2.5. Theoretical constraints 

With so many parameters, it is necessary to restrict the pa- 
rameter space theoretically before comparing the model to 
data. The theoretical constraints arise from multiple con- 
siderations: the (00) Einstein equation, requiring positive 
energy densities, and avoiding instabilities in scalar and 
tensorial perturbations. 

2.5.1. The (00) Einstein equation and eg 

As we use only the (ij) Einstein equation and the tt equa- 
tion of motion to compute the dynamics of the Universe 
(equations |S] and O , we can use the (00) Einstein equation 
as a constraint on the model parameters: 

H 2 = ^ + ^ + £f. i? 2-2_ 2 - 3i j4-3 + 15_ 4i?6 _ 4 _ 7 _ 558 _ 5 ^ 

a d a 4 6 2 

(28) 

More precisely, we use this constraint both at z — to fix 
one of our parameters and, at other redshifts, to check the 
reliability of our numerical computations (see Sect. 13.11) . 
The parameter we choose to fix at z = is: 



with 



1 



c 5 = -(-i + n° m + n° r 



2c, 



15 



-c 4 J 



(29) 



!0 , £2 

6 ' 2 

We choose to fix C5 based on the other parameters because 
allowing it to float introduces significant numerical difficul- 
ties when solving equations |H] and |TJ1 since it represents the 
weight of the most non-linear term in these equations. As 
f2° is fixed given h, our parameter space has been reduced 
to Q^, h, c 2 , c 3 and c 4 . 

2.5.2. Positive energy density 

We require that the energy density of the Galileon field be 
positive from z = to z = 10 7 (see 13 .2. 21 and Appendix B). 
At every redshift in this range, this constraint amounts to: 



^H 2 x 2 -6c 3 H 4 x 3 + ^c i H 6 x 4 -2lc 5 H s x 5 > 0. 

( 3 °) 

This constraint is not really necessary for generic scalar 
field models. But as we will see in the following, it has no 
impact on our analysis as the other theoretical conditions 
described below are stronger. 



2.5.3. Scalar perturbations 

As suggested bv lApplebv k, Linderl (|2012l ). outside the qua- 
sistatic approximation, the propagation equation for Sy 
leads to two conditions, which we again check from z = 
to z = 10 7 to ensure the viability of the linearly perturbed 
model: 

1. a no-ghost condition, which requires a positive energy 
for the perturbation: 



K2 



3k 5 „ 
-— < 

2 K4 



(31) 



k 2 = ~ + 6c 3 H 2 x - 21c±H 4 x 2 + 30c 5 H 6 x 3 . (33) 



2.5.4. Tensorial perturbations 

In this work, we als o add two conditions derived by 
iDe Felice fc Tsuiikawal f|201 lh for the propagation of tensor 
perturbations. Considering a traceless and divergence-free 
perturbation Sgij = a 2 hij, these authors obtain identical 
perturbed actions at second order for each of the two po- 
larisation modes hm and h™. For h ffi : 



SS. 



(2) 



(34) 



with Qt and ct as defined below. From that equation, we 
extract two conditions in our parametrization that have to 
be satisfied (again from z — to z = 10 7 ): 



1. a no- ghost condition 



1 3 

2~t 



4=4 



M 2 = 7; - -ciH x 



c 5 H 5 x 5 > 



(35) 



2. a Laplace stability condition 



\ + \ciH A x A + ^c 5 H b x A {H'x + Hx') 



\c 4 H A x A + §c 5 H 5 x 5 



> 0. (36) 



All the above conditions allow us to reduce our parame- 
ter space significantly. The Galileon model c ontains degen- 
eracie s between the c,*s as pointed out in e.g. IBarreia et al.l 
(2012). The above theoretical constraints and our new 
parametrization allow us to break degeneracies between the 
Ci parameters that would make it difficult to converge to a 
unique best fit with current cosmological observations. As 
an example the tensorial theoretical conditions lead to a 
significant reduction of the parameter space (see dark dot- 
ted regions in Fig. so that closed probability contours 
are obtained. 



3. Likelihood analysis method and observables 

In the following, we define a "scenario" to be a specific 
realisation of the cosmological equations for a given set of 
parameters {fi^, C2, C3, c 4 } . 

To perform t he lik elihood analysis, the method used 
in IConlev et all (|2011l ) for the analysis of SNLS dat£0 
was adapted to the Galileon model. For each cosmologi- 
cal probe, a likelihood surface C is derived by computing 
the x 2 for each visited scenario: £(fi^, C2, C3, c 4 ) cx e~ x / 2 . 
The way h is treated is described in Sect. 13.2.21 Then we 
report the mean value of the marginalized parameters as 
the fit values of 0° and the &s. 



http:/ /casa.colorado.edu/ aaconley/Software.html 
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3.1. Numerical computation method 

In order to compute numerical solutions to equations [5] and 
[9] we use a 4* h order Runge-Kutta method to compute H(z) 
and x(z) iteratively starting from the current epoch, where 
the initial conditions for H and x are specified (see l2.2l) , and 
propagating backwards in time to higher z. We use a suffi- 
ciently small step size in z to avoid numerical divergences. 
This is challenging because of the significant non-linearities 
in our equations. To determine the step size, we therefore 
require that equation |2"B1 normalized by H 2 , be satisfied at 
better than 10~ 5 for each step. 

At each step of the computation, we also check that all 
the previously discussed theoretical conditions are satisfied 
(equations [30] [3TJ [32] [35] and [36]) . Cosmological scenarios 
that fail any of these conditions are rejected and their likeli- 
hood set to zero. The result of these requirements is shown 
e.g. in Fig.[2]as dark dotted regions. Eauationl30lconcerns a 
negligible number of Galileon scenarios, but the four other 
constraints lead to a significant reduction of the parameter 
space. 

3.2. Data 

Here we describe the cosmological observations that we 
make use of in our analysis. Special care was taken to 
choose data that do not depend on additional cosmolog- 
ical assumptions. 

3.2.1. Type la supernovae 

The SN la data sa mple used in this w ork is the SNLS3 
sample described in lConlev et al.l (|2011l) . It consists of 472 
well measured supernovae from the SNLS, SDSS, HST, and 
a variety of low-z surveys. 

A type la supernova with intrinsic stretch s and color 
C has a rest-frame B-band apparent magnitude raj which 
can be modeled as follows: 

m™ od = 5log w V L (zhei,zcMB,cosmo)-a(s-l)+0.C+MB 

(37) 

T>l is the Hubble-constant free luminosity distance, which 
in a flat Universe is given by: 



£> L(zheh ZCMB, COSmo) = (1 + Z hel ) 



dz 



H(z, cosmo) 

Zhei and zcmb are the SN la redshift in the heliocentric 
and CMB rest frames respectively, "cosmo" represents the 
cosmological parameters of the model, a and (3 are pa- 
rameters describing the light-curve width-luminosity and 
colour-luminosity relationships for SNe la. Mb is defined 
as Mb = Mb + 51og 10 c/H + 25, where Mb is the rest- 
frame absolute magnitude of a fiducial (s = 1, C = 0) SN la 
in the B-band, and c/H is expressed in Mpc. a, j3 and Mb 
are nuisance parameters which ar e fit simultaneously with 
the cosmological par ameters. As in lConlev et al.l f|2011h and 
ISullivan et atl(|2011lV we allow for different M b in galaxies 
with the host galaxy stellar mass below and above 10 10 M Q 
to account for relations between SN la brightness and host 
properties that are not corrected for via the standard s and 
C relations. When computing type la supernova distance lu- 
minosities in Sect. [H we neglect the radiation component 
in H(z), since all measurements are restricted to redshifts 



below 1.4 where the effects of radiation density are negligi- 
ble. 

Systematic uncertainties must be treated carefully when 
using SN la data, due to their dependence on a and /? and 
due to covariances between diffe rent supernovae. In this 
paper we follow the treatment of IConlev et al.l (]2011| ) and 
ISullivan et aD (1201 lb . 

3.2.2. Cosmological Microwave Background 

The CMB is a powerful probe to constraint the expansion 
history of the Universe as it gives high redshift cosmolog- 
ical observables. The power spectrum provides a lot of in- 
formation on the content of the Universe and the relations 
between the different fluids, as long as we are able to model 
the thermodynamics of these fluids before recombination. 
The Galileon model does not modify the standard baryon- 
photon flux physics as long as the Galileon field docs not 
couple directly to matter, as is assumed in this work. Thus, 
the usual formulae and predictions used in the standard 
analysis of the CMB power spectrum remain valid. 

The positions of the acoustic peaks ca n be quantified 
by th ree o bservables: Mn, R, z*) (see e.g. iKomatsu et al.l 
(|201lD and IKomatsu et ah! ()2009D ) where l a is the acous- 
tic scale related to the comoving sound speed horizon, R is 
the shift parameter related to the distance between us and 
the last scattering surface, and z* is the redshift of the last 
scattering surface. These quantities are derived from the 
angular diameter distance, which in a flat space is given 
by: 

^ / n c 1 f z dz' 
D a{z) = - e - t — / (39) 



(40) 



H l + zJ H{z>) 
and from the comoving sound speed horizon: 

r s (z) = jr / da - 

H J a 2 H(a) 



c s is the usual normalized sound speed in the baryon- 
photon fluid before recombination: 



1 



l + 3(3r*£/4fi°)a 



(41) 



where is the baryon energy density parameter today. 
With the above definitions, the acoustic scale l a is given 

by: 

l. = (l + *.)^%) (42) 
and the shift parameter R by: 

fl=s ( i +#A(2 , )= ^r|_, ( 43) 



is given by the fitting formula of lHu fc Sugivarrial (|l996h : 

z* = 1048 [l + 0.00124(^ 2 )- - 738 ] [l + gi (n° n h 2 ) g2 ] 

(44) 



.91 



.92 = 



0.0783(»gfr 2 )-°- 
1 + 39.5(f^ 2 ) - 763 

0.560 
l + 21.1(ft^ 2 ) 1 - 81 



(45) 
(46) 
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According to iHu fc Sueivamal (|l996t) , formula HH is valid 
for a wide range of ^Jt 2 and f2°/i 2 . 

In order to compare these observables to the 7-year 
WMAP data (WMAP7), we follow the numerical recipe 
given m iKomatsu et~aT1 J2009). The key point of this recipe 
is that, for each cosmological scenario, Xcmb mus t be min- 
imized over h and £^/i 2 , which appear in equation 1441 and 
in the computation of H(z) through (see equation [l"4"f . 

An important feature to note is that we have to solve 
equations [5] and [5] from a = 1 to a = to compute the 
CMB observables. Numerically, however, we can't reach a = 
[z = oo) due to numerical divergences. To avoid such 
numerical complications, we carry out these computations 
up to a = 10~ 7 and then linearly extrapolate the value of 
the integral to a = (for more details on the reliability of 
this approximation see Appendix B). Thus, the theoretical 
constraints of 12. 51 are checked from a = 1 to a = 10~ 7 . 

Finally, as CMB observables depend explicitly on H , 
we impose a G aussian prior on its value, h = 0.737 ± 0.024 
as measured bv lRiess et al.l ((201 lh from low redshift SNe la 
and Cepheid variables. 

The WMAP 7 recommended best fit values of the CMB 
observables are: 



/ (L) 

(Vcmb) = (R) 



302.09 ±0.76 ' 
1.725 ±0.018 
1091.3 ±0.91 



with the corresponding inverse covariance matrix: 



J CMB 



( 2.305 29.698 
= 29.698 6825.270 
V -1.333 -113.180 



-1.333 1 
113.180 
3.414 



(47) 



(48) 



from Komatsu et al.l (|2011l) . As pointed out by 
iNesseris et al.l i 20101 ). the unco upled Galileon mode l 
fulfils the assumptions required in IKomatsu et al.l ([2009) 
to use these distance priors, namely a FLRW Universe with 
the standard number of neutrinos and a dark energy back- 
ground with negligible interactions with the primordial 
Universe. Once the observables {l a , R, z*} are computed in 
a cosmological scenario, we build the difference vector: 



A Vcmb = R 
\ z~ 



-(V. 



CMB) 



and compute the CMB contribution to the total x 2 : 

(fr- 0.738) 2 



XCMB+H„ - AV CMB C CMB AV C MB- 



0.024 2 



(49) 



(50) 



3.2.3. Baryonic Acoustic Oscillations 

BAO distances provide information on the imprint of the 
comoving sound horizon after recombination on the dis- 
tribution of galaxies. The BAO observable is defined as 
y s (z) = r s {zd) I 'Dv(z) where r s is the comoving sound hori- 
zon at the baryon drag epoch redshift Zd, and Dy(z) is the 
effective distance (Eiscnste in et al.l . 120051 ) given by: 



D v (z) 



"I 1/3 



{l + zfD\(z) 



H(z) 



(51) 



Zd is computed using the lEisenstein fc Hul (|l998l ) fitting 
formula: 



Zd 



i29i(r;°„fr 2 ) - 251 

1 + 0.659(0°, ^2)0.828 



(52) 



h = 0.3l3(n° m h 2 y OA19 [1 + 0.607(f> m /i 2 ) - 674 ] 
b 2 = 0.238(n° m h 2 ) - 223 . 



(53) 
(54) 



This formula remains valid in presence of a Galileon field 
not coupled to matter. 

As for the CMB, BAO distances depend on h and so 
we follow the same recipe as previously mentioned to com - 
pute them, including the H prior from lRiess et al.l (|201ll ). 
We also make the same approximation to compute r s as 
for the CMB. The minimization over h and il®h 2 is done 
independently for CMB and BAO when their individual 
constraints are derived and simultaneously when combined 
constraints are computed. 

We use the dataset of distances derived from galaxy sur- 
veys as published in the SD SS-III BOSS cosmological a nal- 
ysis (|Anderson et al.l ([2012D and I Sanchez et all (|2012t )) to 
avoid redshift overlaps in the measurements (see Table [T]). 

For a cosmological constraint derived from BAO dis- 
tances alone, the BAO contribution to the total x 2 is given 
by: 



Xbao+Hq 



E 



(y s (z)-y? es (z)) 2 , (fe- 0.738) 2 



0.024 2 

{9P b h 2 - 0.02249) 2 
0.00057 2 



(55) 



where we have added a Gaussian prior on Q®h 2 when deal- 
ing with this probe alone. 

When BAO and CMB probes are combined, we compute 
their contributions to the x 2 simultaneously in order to 
not over-count the Hubble constant prior. Therefore, the 
combined contribution is: 



XCMB+BAO+H c 



- AV CA/B C CMB AV CAfS 

-yT es (z)) 2 (Zi- 0.738) 2 



(56) 



0.024 2 



3.2.4. Growth rate of structures 



The cosmological growth of structure is a critical test of 
the G alileon model, as noted by many authors (see iLinderl 
(2005) for example). It is a very discriminant constraint for 
distinguishing dark energy and modified gravity models. 
Many models can mimic ACDM behaviour for the expan- 
sion history of the Universe, but all modify gravity and 
structure formation in a different manner. 

In linear perturbation theory, the growth of a matter 
perturbation S m = 5p m /p m is governed by the equation: 



2HS„ 



4nG NPm 6 m = 0. (57) 

But a s argued in ILinderl (|2005[ ) and used in IKomatsu et al.l 
(|2009l) . it is better to study the growth evolution with the 
function g(a) = D(a)/a = 6 m (a) / (aS m (l)). In the Galileon 

case, the Newton's constant is replaced by (a) as given 
in equation[57] The g{a) is obtained by solving the following 
second order differential equation: 



d 2 g 1 / a dH 
_| 5 H = 

da 2 a \ H da 
1 



dg 
da 

a dH 3 G 



s 3+^— -- 



(</>) 

eff 



H da 2 G N a 3 H 2 



0. (58) 
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Table 1. BAO measurements used in this work. 



z 




Survey 


Reference 


0.106 


0.336 ±0.015 


6dFGS 


Beutler et al. (2011) 


0.35 


0.1126 ±0.0022 


SDSS LRG 


Padmanabhan et al. ('2012) 


0.57 


0.0732 ± 0.0012 


BOSS CMASS 


Anderson et al. (2012) 



Table 2. Growth data used in this work. 



z 


fo%(z) 


F(z) 


r 


Survey 


Reference 


0.067 


0.423 ± 0.055 






6dFGRS (a) 


Beutler et al. (2012) 


0.17 


0.51 ±0.06 






2dFGRS (a) 


Percival et al. (2004) 


0.22 


0.53 ±0.14 


0.28 ± 0.04 


0.83 


WiggleZ 


Blake et al. ( 201 lb) 


0.25 


0.351 ±0.058 






SDSS LRG (b) 


Samushia et al. (2012a) 


0.37 


0.460 ± 0.038 






SDSS LRG (b) 


Samushia et al. (2012a) 


0.41 


0.40 ±0.13 


0.44 ± 0.07 


0.94 


WiggleZ 


Blake et al. (2011b) 


0.57 


0.430 ± 0.067 


0.677 ± 0.042 


0.871 


BOSS CMASS 


Reid et al. (2012) 


0.6 


0.37 ±0.08 


0.68 ±0.06 


0.89 


WiggleZ 


Blake et al. (2011b) 


0.78 


0.49 ±0.12 


0.49 ±0.12 


0.84 


WiggleZ 


Blake et al. (2011b) 



Notes, r is the cross-correlation in (F, fas)- (a) Alcock-Paczynski effect is negligible at low redshift. (b) Values of fag are corrected 
for the Alcock-Paczynski effect but no F(z) values are provided. 



Galileon presence. This equation holds if D(a) has no scale 
dependence which is the case in both models in the linear 
regime. Equation l59l takes into account the different growth 
histories since recombination in the two models. 

However, standalone fas(a) measurements extracted 
from observed matter power spectra usually make use of a 
fiducial cosmology which assumes General Relativity. This 
hypothesis is no longer neces sary when taking into acco unt 
the Alcock-Paczynski effect (|Alcock fc Paczvnskil Il979t ) in 
the power spectrum analysis. This results in joint mea- 
surements of f<Ts(a) and the Alcock-Paczynski parameter 
F (a) = c~ 1 DA(fl)H(a)/a which are to be p referred when 
constr aini ng modified gravi t y mode ls (see e.g. lBeutler et all 
(2012) and lSamushia et al.l (j2012bT )V Note that equations^] 
and[5]are all we need to predict F(a) in the Galileon model. 

The measurements of fas(z) and F{z) used in this work 
are summarized in Table [5] To compare these with our 
model, we first solve equations [5] and |H] from a = 1 to ainitiai 
to obtain values of H(a), F(a) and G^g (a)/Gjv, and then 
solve equation [55] from ainitiai to a = 1 which provides us 
with fas{z) predictions. 

As F{z) and fas(z) measurements are correlated, a 
covariance matrix CqoS is built using data presented in 
Tabled Moreover our /erg prediction relies on the WMAP7 
measurement of cr 8 (a = 1) (equation [60]), so the WMAP7 
experimental uncertainty is also propagated to the diago- 
nal and off-diagonal terms of CgoS- Then a vector VgoS 



A natural choice for t he initial conditions is (/(ainitiai) = 1 
and dg/da |a initial = ([Komatsu et al.1 [20091) . where owia! 
is 0.001 ~ 1/(1 + z*). We checked that our results do not 
depend on this choice as long as ainitiai is taken between 
10~ 2 and 10~ 5 . 

Measurements of the rate of growth of cosmic struc- 
tures from redshift space distortions can be expressed in 
terms of /(a) = dlnD(a)/dlna or /crg(a), where erg is 
the normalization of the matter power spectrum. fas(a) 
is known to be less sensitive to the overall normalization 
of the power spectrum model used to derive the measure- 
ments (jSong fc Percivalll2009D . So this is the observable we 
choose in this work. To make predictions for /crs(a) in our 
analysis, we solve equation [35] to obtain g(a), from which 
we deduce f( a) and D(a), and we co mpute os(a) in the 
following way ([Samushia et all l2012al) : 

D(a) 

08 (a) = cr 8 (ainitiai) 7^7 r (59) 

^(ainitiai) 

where 

riACDM / \ 

0-8 (flinitial) - O-g [a - 1 ) TjACDM?^— T) ( 60 ) 

and <rf MAP7 (a = 1) = Q.811±g;gl? is the present value 
of the CMB power sp ectrum normalization published by 
iKomatsu et al.l ([20 111) in the framework of the ACDM 
model. Equation [50] states that the normalization of the 
CMB power spectrum at decoupling is the same in the 
ACDM and Galileon models, which is consistent with our 
assumption that the CMB physics is not modified by the 
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containing all predictions at each z% is built: 



V GoS = 



\ 



J 



(61) 



The contribution of the growth rate of structure to the total 

x 2 



2 

XGoS 



(62) 



with AV GoS = V G oS - (V Go s) where (V GoS ) contains the 
measurements of Table [2J 

Note that equation Q3] requires a value for f2°, and hence 
in principle this equation should be simultaneously solved 
with the BAO and CMB constraints using the same prior 
on Hq. However, we find that doing so has essentially no 
effect on our \ 2 ■ So, to accelerate the computation, we set 
here h to t h e valu e derived from the Ho measurements of 
iRiess et all (|2011l ). 

4. Results 

In the following we present the results of the experimental 
constraints on the Galileon model derived from the cosmo- 
logical probes. 

4.1. SN constraints 

Results from SN la data are presented in Fig.[5]and TableEl 

4.1.1. SN results 

Despite the large number of free parameters in the model, 
we obtained closed probability contours in any two- 
dimensional projection of the parameter space. We observe 
strong correlations between the 5jS, especially between C2 
and C3. 

We note that the best fit value for Q^n ~ 0-27 is com- 
patible with the current constraints obtained in the ACDM 
or FWCDM models. The CjS are found to be globally 
of order ss —1. From the best fit values of the parame- 
ters, we derive the value of C5 using equation [2H1 and find 
c 5 = —0. 349+^555, including systematic uncertainties. 

In the following we discuss the impact of fixing the nui- 
sance parameters a and (3 and the effect of systematics on 
the best fit values. 

4.1.2. Impact of nuisance parameters 

When marginalizing over the cosmological parameters, the 
best fit values of the SN nuisance parameters a, /3, Ai B and 
A4 B in the Galileon context are identical to those published 
for the ACDM model, as shown in Fig.ffland TableU This is 
a truly important point to note. It means that the modelling 
of the SN la physics contained in these nuisance parameters 
is adequate for these two cosmological models despite their 
differences. 

In principle, the correct method to use when analyzing 
SN la data is to scan and marginalize over the nuisance 
parameters. However, once the best fit values of a and /? 
are known, keeping them fixed to their best fit values in any 



3.6 - 



3.4 - 



3.2 



3.0 - 



2.8 



2.6 - 




1.0 1.2 



1.6 



1.8 2.0 



Fig. 1. Confidence contours for the SN nuisance parameters 
a and /3 when marginalizing over all other parameters of 
the model. Dashed red contours represent 68.3%, 95.4% and 
99.7% probability contours for the ACDM model. Filled 
blue contours are for the Galileon model. Note that they 
are nearly identical, the Galileon one is just 2.8% wider 
which is likely due to larger steps in a and /3. See Table [3] 
for numerical values. 



study using the same SN sample has a negligible impact on 
our results (see Table [3]) . In the Galileon case, the contour 
areas decrease only by 0.7% and have the same shape as 
in Fig. [2j So, for further studies with the SNLS3 sample in 
the ACDM or Galileon models, our analysis demonstrates 
that it is reasonable to keep the nuisance parameters fixed 
to their published values in the SNLS papers. 

4.1.3. Impact of systematic uncertainties 

From the results in Tabled we note that the identified sys- 
tematic uncertainties shift the best fit values of the Galileon 
parameters by less than their statistical uncertainties. With 
systematics included the area of the inner contours in- 
creases by about 53%. This is less than what is observed in 
fits to the ACD M or Flat wCDM m odels (103% and 80% 
respectively, see IConlev et al] (|2011h V 

4.2. Combined CMB, BAO and H constraints 

The results using CMB, BAO and Hq data are presented 
in Fig. [3] and Table g] 

The combined WMAP7+BAO+H0 data provide a very 
powerful constraint on fi^j, but no tighter constraints on 



the & than SNe la alone, ilz 



0.272 



+0.014 
-0.009 



is, as for the 



SNLS3 sample, close to the current best estimates for this 
parameter in the standard cosmologies, but this time with 
very sharp error bars competitive with the most recent 
studies on other cosmological models. However, the Cj best 
fit values are similar to those predicted with the SNLS3 
sample. 

To use the WMAP7+BAO+H0 data, h and n° b h 2 
have to be minimized for each explored Galileon scenario. 
Minimized values of these parameters are collected in the 
histograms of Fig. [3] for the subset of the scenarios that 
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Fig. 2. Experimental constraints on the Galileon model from SNLS3 data alone. To represent the four-dimensional likeli- 
hood £(f2^ i; C2, C3, C4), six two-dimensional contours for each pair of the Galileon model parameters are presented, after 
marginalizing over JvO^, cv, /?, and the remaining Galileon parameters. The filled dark, medium and light blue 

contours enclose 68.3, 95.4 and 99.7% of the probability, respectively. The contours include statistical and all identified 
systematic uncertainties. The dark dotted regions correspond to scenarios rejected by theoretical constraints, as described 
in the text. Labels in these regions indicate the main cause for excluding the scenarios. 



Table 3. Cosmological constraints on the Galileon model from the SNLS3 sample 



Method 


n° 


C2 


C3 




f'4 


a 


P 


Mi, 


Ml 


x 2 


Stat+sys+a/3 


n 97 o+0.057 
U.ZI O_ 42 




- 1.779 j 


-1.073 
-1.416 


-0.587l°j 4 « 


1 49S+ ' 121 


q o/?q+0.121 
o - iuJ -0 .103 


23.997 


23.950 


415.4 


Stat+sys 


n 97 o+0.054 
<J.ZIO_ 042 


_c; o4n +1 ' 880 
J.Z'iU_2.go2 


-1.78lj 


-1.071 
-1.426 


-0.588l°:£| 


1.428 


3.263 


23.997 


23.950 


420.1 


Stat only 


294+ 0045 

U.Z»4_ 039 


*-' DO -2.921 


-1.586j 


-0.987 
-1.474 


-0.541«;f 3 2 8 


1.451 


3.165 


24.022 


23.951 


441.8 



Notes. Results are computed using either statistical and systematic uncertainties combined or statistical uncertainties only. In 
the first line, we marginalize over a and f3, whereas in the last two lines, a and (3 are kept fixed to t heir marginalized values . No 
errors are given on Mg and M% as they are analytically marginalized over fsee lConlev etHI IpoTlT )). 
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0.2 0.3 



Fig. 3. Experimental constraints on the Galileon model from WMAP7+BAO+H0 data. To represent the four-dimensional 
likelihood £(f2^, c~2, C3, C4), six two-dimensional contours for each pair of the Galileon model parameters are presented 
after marginalizing over the left over Galileon parameters. The filled dark, medium and light green contours enclose 
68.3, 95.4 and 99.7% of the probability, respectively. Dark dotted regions correspond to scenarios rejected by theoretical 
constraints. 



fulfilled the theoretical constraints. Values for the best fit 
scenarios are reported in Table 2] For the Galileon model, 
the h distribution has a mean of 0.65 with a dispersion of 
0.06, compatible with th e Hn prior. The c onstraint on h 
is slightly lower than the lRiess et al.l (|2011l ) value, but the 
same behaviour is obtained for the ACDM model using the 
same program and data. The central value for the Q®h 2 
distribution is fully compatible with the WMAP7 value, in 
both the Galileon model and the ACDM model. However, 
in the Galileon model, values below 0.22 are much more 
disfavoured. 

For completeness, we present in Fig. [5] examples of re- 
sults obtained from the WMAP7+H0 and BAO+H0 probes 
separately. Both plots have been obtained with a minimiza- 
tion on h and Sl^/z 2 , but a Gaussian prior on fl®h 2 was 
added for the BAO (see Equation [55]). We use the WMAP7 
constraint for that prior since Fig.2]shows that the Galileon 
model is consistent with it. 



4.3. Growth of structure constraints 

Results using growth data are presented in Fig. [5] and 
Table HI and further commented in Sect. [5] 

Growth data and cosmological distances provide consis- 
tent values for the c$s. The best fit value from growth 
data, Sl^ w 0.20, is below that from the other probes but 
still compatible at the 1.5er level. This is the main difference 
between the two types of probes. 

However, the use of growth data in cosmology deserves 
some comments. In our work as in many others, differ- 
ent assumptions about the importance of non-linearities in 
structure formation are made in the theoretical predictions 
and in the experimental extraction of growth data from the 
measured matter power spectrum. 

As noted in Sect. 12.41 our theoretical predictions are 
derived in the line ar regime, and using a quasistatic ap- 
proximation. While iBarreia et~ al, (2012) confirm that the 
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0.02 0.021 0.022 0.023 0.024 0.02S 0.026 0.027 0.028 0.029 



n n£h 2 

Fig. 4. Minimized values of h and £l®h 2 for a large sub- 
set of tested scenarios, in ACDM (red dashed histogram) 
and in the Galileon cosmology (blue filled histogram). 
Dashed black band s represent the m easurements of Hp from 
IRiess et all (|2011ft and n°h 2 from iKomatsu et all (j2011ft . 
Only scenarios with % 2 < 200 enter these histograms to 
deal only with pertinent scenarios. Note that both models 
give values of h and fl®h 2 in agreement with the measure- 
ments. 



latter is valid in the Galileon model, using only the linear 
regime is restrictive. As an example, this may be the ori- 
gin of the divergences in G^(z )/Gn that appear i n som e 
Galileon scenarios as noted by lApplebv fc Linden (|2012ft . 
Going beyond linear perturbation theory may change our 
predictions and thus could modify the result of our analysis. 

In order to estimate the previous effect, we can try to 
identify at which scale non-linearities start to matter and to 
check whether this value is outside the range of scales taken 
into account in the growth of structure measurements. As 
an example, WiggleZ measurements of fas are derived us- 
ing a non-linear growth of structure model ([Jennings et al.l . 
I2011h and encompass all scales k < 0.3 hMpc . In this 
model, the frontier between the linear and the non-linear 
regimes is k w 0.03 /iMpc -1 . Other measurements in 
Table [2] include scales up to k ss 0.2-0.4 /iMpc -1 too. 
On the other hand, there is no prediction in the Galileon 
model which goes beyond the linear regime. However, esti- 
mates of the scale at which non-linear effects appear exist 
in similar modified gravity models. Numerical simulations 
of the Chameleon screening effect for / (R) theories show 
that non-linearit y effects can be important at scales k sa 
0.05 hMpc- 1 (seelBrax et alj (|2012ft . l Jennings et al.l ()2012ft 



and iLi et all (12012ft '). Other simulations of the Vainshtein 
effect in the DGP model show that significant differences 
between the linea r and non-linea r regimes appear for scales 
k > 0.2 feMpc- 1 (|Schmidtl . [2009ft . 

So the lack of non-linear effects in our perturbation 
equations, and hence in our predictions for /erg in the 
Galileon model, is likely to have a significant impact on the 
constraints we derived from growth measurements, since 
the latter accounted partially for non-linear effects. 

4.4. Full combined constraints 

Results from all data are presented in Fig. [7]Table 0] 
presents the best fit values for the Galileon model parame- 





Fig. 5. Experimental constraints on the Galileon parame- 
ters and c 2 from WMAP7+H0 data (top panel) and 
from BAO+H0 data (bottom panel). The four-dimensional 
likelihood £(^1^, C2, S3, C4) has been marginalized over C3 
and C4. The filled dark, medium and light green contours 
enclose 68.3, 95.4 and 99.7% of the probability, respectively. 
Dark dotted regions correspond to scenarios rejected by 
theoretical constraints. A Gaussian prior on 51° h 2 based on 
the WMAP7 value has been added to the BAO+H0 fit. 



ters. The derived C5 value is: 



-best fit 



-0 ^78+ - 120 

u -°' 8 -0.219- 



(63) 



Note that negative values are preferred for the CjS at the la 
level. M oreover, the Galileo n h best fit values is compatible 
with the IRiess et~aTI (|2011f ) measurement. 

We carried out an a posteriori check to identify which 
scenarios present a significant amount of early dark energy. 
At decoupling, fi^z*) > 10%Q r (z*) only for viable scenar- 
ios with Q^, > 0.3, c 2 > -4, c 3 > -1 and c 4 > 0. Note 
that this check can be done after comparing theory with 
data as only data can provide values for h and il®h 2 . So, 
for Galileon scenarios with 51^ < 0-3, which is the region 
favoured by data, we find no significant early dark energy. 
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0.5 



0.0 0.1 0.2 0.3 0.4 0.5 -12 -10 -8 -6 -4 -2 

n m c2 

Fig. 6. Experimental constraints on the Galileon model from growth data (red) and from SNLS3+WMAP7+BAO+H0 
combined constraints (dashed). The filled dark, medium and light coloured contours enclose 68.3, 95.4 and 99.7% of the 
probability, respectively. Dark dotted regions correspond to scenarios rejected by theoretical constraints. 



4.5. Analysis of the best fit scenario 



4.5.1. Cosmic evolution 



What does the best fit scenario (derived from all of the data; 
the last line of table 2]) look like? As p^ can be defined from 
the (00) Einstein equation, a Galileon pressure can be 
defined from the (ij) Einstein equation: 



Fx 



H 2 M 2 2 



c 4 



^-H 2 x 2 + 2c 3 H 3 x 2 (Hx)' 



9 



-H b x* + l2H b x 6 x' + 15H°x A H' 



+ 3c 5 H 7 x 4 (5Hx' + 7H'x + 2Hx) . (64) 

Combining p w and P„-, an equation of state parameter 
w 7r {z) = Pn{z) I P-k{z) can be built for the Galileon "fluid". 
We can also construct an equation for ^{z) using p v {z) — 
D, 7r (z)HQMp/(3H 2 (z)). The evolution of w % (z), Q„(z) and 

(z) /G N for the Galileon best fit scenario are shown in 
Figs. [S] and |U 



The left panel of Fig. [S] shows that, for the best fit scenario, 
radiation, matter and dark energy (here the Galileon) dom- 
inate alternatively during the history of the Universe, as in 
any standard cosmological model. These three epochs are 
also visible in the evolution of w(z). Moreover, the best fit 
scenario evolves in the future towards the de Sitter solu- 
tion w = — 1, which is a n attr actor of the Galileon model 
(|De Felice fc TsuiikawaL l2010D . In the region < z < 1, 
where SNe put a tight constraint on dark energy, w(z) de- 
viates significantly from -1, its ACDM value. Note that in 
the fit with SNe alone, the deviation is less pronounced, 
with an average value of -1.09 in < z < 1 which is com- 
patible with the fitted v alue of w in constant w dark energy 
models, as published in lConlev et al.l (|201ll) . 

During matter domination, dark energy contributes 
about 0.4% to the mass-energy budget at z = 10. For 
comparison, in a standard ACDM model dark energy con- 
tributes only 0.2% at this redshift (assuming a flat ACDM 
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Fig. 7. Combined constraints on the Galileon model from SNLS3, WMAP7+BAO+H0 and growth data. The filled dark, 
medium and light yellow contours enclose 68.3, 95.4 and 99.7% of the probability, respectively. Dark dotted regions 
correspond to scenarios rejected by theoretical constraints. 
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Fig. 8. Evolution of the Sli(z) (left) and of w(z) (right, solid curve) for the best fit Galileon model from all data (last 
row of Tabled]). As a comparison, the dashed orange line gives w(z) for the best fit scenario from SN data alone. 
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Table 4. Galileon model best fit values from different data samples 



Probe 




C2 


C3 


C4 


/) 


n° b h 2 


x 2 


SNLS3 


n 070+0. 054 
UJ '°-0.042 


-5.240±J : »° 


-1-'81_1.426 


-0.588±°j 4 6 8 


. 


_ 


420.1 


Growth 


n 2nn+ 047 

u.zuu_ 044 


-5.430i;;«S 


1 7C7+0.365 
— J -'°'-1.251 


-0.6351°,;^ 






19.83 


BAO+WMAP7+H0 


n 070+0.014 

U.z, 1 ^_ .009 


~~ 5.591 1 2 .655 


-1 926+ 1 008 
i.»zu_ 1 407 


-0-619l°;f 3 8 5 


0.713 


0.0224 


2.14 


SNLS3+BAO+WMAP7+H0 


n 070+O.O14 
U.z, 1 ^_ .008 


O.ODJ_ 2 gg 4 


-1 Q17+ 1001 


-0.619l°;f 3 8 3 


0.713 


0.0224 


423.1 


SNLS3+BAO+WMAP7+H0+Growth 


n 071 +0.013 

K, - L 1 i-o.oos 


4 nro + 0.518 


1 KQ7+0.203 
X.oai -0.726 


n 77-1 +0.098 
<J- 1 1 r_o.o6l 


0.735 


0.0220 


450.4 



Notes. SNLS3 with systematics included, a and /3 fixed to their marginalized value, h and Q®h 2 have been minimized so no error 
bars are provided. 



Table 5. ACDM best fit values from different data samples 



Probe 






h 


n°h 2 


x 2 


SNLS3 


fl 1 78+0-100 


0-6641°:™ 






419.7 


Growth 


0.295ir 3 3 I 


646+ 067 

u -°^°-0.072 






8.2 


BAO+WMAP7+H0 


0.2881°;^ 


n 71 o+0. 016 
u. 1 rO-o.014 


0.691 


0.0225 


5.6 


SNLS3+BAO+WMAP7+H0 


n oQQ+0.013 
"• zoo -o.oio 


71Q+0.016 

U. 1 J-»- .oi3 


0.692 


0.0225 


427.8 


SNLS3+BAO+WMAP7+H0+Growth 


n 077+0.011 

u - z ' ' -0.009 


fl 7oc;+0.015 
U. <^0_o .012 


0.698 


0.0225 


440.2 



Notes. SNLS3 with systematics included, a and /3 fixed to their marginalized value, h and f2°/i 2 have been minimized so no error 
bars are provided. 



model with fij^ = 0.27). In the same way, dark energy 
contributes 0.04% at z„ in the Galileon best fit scenario, 
whereas for ACDM J7a = 10~ 9 at z*. In our best fit Galileon 
scenario, dark energy is more present throughout the his- 
tory of the Universe than in the ACDM model, but is still 
negligible during the matter and radiation eras. 

Figure |9] shows the evolution of G^{z)/Gn for the 
best fit scenario and for the growth data best fit scenario. 
Both curves show deviations from 1 at redshifts around 0. 
Particularly, the divergence near the current epoch suggests 
that we should push the Galileon predictions for fcr% be- 
yond the linear regime as already advocated in Sect. 14.31 



4.5.2. Comparison with ACDM 

In Fig. Q3H and Table [5J best fit values for the ACDM pa- 
rameters are presented using the same analysis tools and 
observables. Interestingly, even in the ACDM model there 
is tension between growth data and other probes. The f2J^ 
best fit value is similar in both mode l s, but the h value de- 
parts more from the H iRiess et ail (|2011|) measurement. 
As far as the x 2 s are concerned, SNe la provide a good 
agreement with both models. CMB+BAO+H0 data are 
more compatible with the Galileon model reflecting the 
better agreement on the h minimized value. Yet growth 



U 



100 

l+z 



Fig. 9. Evolution of G ( c f(z)/G N for the best fit scenario 
from growth data only (dashed orange line) and from all 
data (blue solid line). 



of structure data are in better agreement with the ACDM 
model. Finally, due to the poorer fit to growth data in the 
Galileon model, the difference in \ 2 is A^ 2 = 10.2. This in- 
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dicates that the Galileon model is slightly disfavoured with 
respect to the ACDM model, despite having 2 extra free 
parameters. 

As we are comparing two models with different number 
of parameters and complexity, other criteria than compar- 
ing x 2 s can be helpfu l . A re view of selection model critcris 
is provided in iLiddld ((2007) . As our study leads to the full 
computation of the likelihood functions, w e can use pre- 



summation over all 9s: 



cise c ri teria such as the Baye s factor (see Beringer et al.l 



(IMl . ljohn k Narhkari(l2^0a . lKas7fc Raftervl (| 19951) and 



Liddlc (2009) ) or the Devianc e Inf ormation Cr i terion (DIC, 
see lSpiegelhalter et al.l (120021) andlKunz et all (120061) ). AIC 
and BIC criteria used in Nesseris et al.l ( 201C ) are approx- 
imations of the first two using only the maximum of the 
likelihood and not the whole function. Hereafter we restrict 
the discussion to the DIC criterion. 




Fig. 10. Experimental constraints on the ACDM 
model from SNLS3 data (blue), growth data (red), 
BAO+WMAP7+H0 data (green) and all data combined 
(yellow). Black dashed line indicates the flatness condition 
TO + Oa = 1 . 



The DIC criterion is based on the computation of the 
deviance likelihoods Dev(9) = —2\ogp(D \ 9) + C (with 
C a constant not important for DIC evaluation). p(D \ 9) 
is the computed likelihood function C(9) of the model. An 
effective number of parameters po — Dev(9) — Dev(9) is 
derived with 9 the expectation values for 9 and Dev(#) the 
mean deviance likelihood value: 



Dev(0) = -2 J d9p{9 | D)log£(0) 



(65) 



p{9 | D) is the posterior probability density function for a 
vector 9 of parameters of the tested model, knowing the 
data D: 

p(D | 9) x prior(6») 



p(9 | D) 



p(D) 



(66) 



p(D) = J d9 p(D | 9) x prior((9) = J d9C{9) x prior((9) 

(67) 

Note that if the priors are flat, p(9 | D) is just the likeli- 
hood function C{9) normalized to 1. In our case, prior(6*) 
is a flat prior reflecting the theoretically allowed volume 
in the scanned parameter space. We checked that the DIC 
criterion is not sensitive to the exact definition of the prior, 
which makes it a robust tool. 

Then DIC = Dev(0) + 2p D = Dcv(<9) +p D . The model 
with the smallest DIC is favoured by data. In our study, 
we obtain DICGaiiicon — DICacdm = 12.25 > 0. Again, the 
Galileon model is slightly disfavoured by data against the 
ACDM model. The DIC criterion just reflects the A\ 2 and 
does not penalise so much the Galileon model due to its 
higher number of free parameters. 

In the future, provided the tension between growth of 
structure data and distances is not worse after more pre- 
cise measurements of the observables used in this paper are 
included, new observables will be necessary to disentangle 
between the two models. A promising wa y would be to ex- 
ploit e.g. the ISW effect as discussed in iKobavashi et al.l 
(I2010D . 

4.5.3. Comparison with FWCDM 




0.6 



p(D), the probability to obtain the data D, is also called 
the marginal likelihood as it can be computed using the 



Fig. 11. Experimental constraints on the FWCDM 
model from SNLS3 data (blue), growth data (red), 
BAO+WMAP7+H0 data (green) and all data combined 
(yellow) . 



For consistency with our assumption about flatness, 
we present also a comparison with the effective FWCDM 
model, a model with constant dark energy equation of state 
parameter w in a flat Universe (see Table [6] and Fig. [TTj) . 
The data set points toward a value of w below -1, which is 
consistent with the Galileon best fit scenario (see Fig. [5]). 

However, the difference in x 2 is the same as for the 
ACDM model, A\ 2 = 10.2, and the DIC criterion gives 
DIC G aiiicon - DICfwcdm = 12.16 > 0. Here again, the 
Galileon model is not significantly disfavoured. 
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Table 6. FWCDM best fit values from different data samples 



Probe 




w 


h 


Q° b h 2 


x 2 


SNLS3 


0.1831°;™ 


-0 Q1+ 017 


. 


. 


419.6 


Growth 


0.294t° ;°» 


n o 7 +0. 09 






7.9 


BAO+WMAP7+H0 


n 9 77 +0.017 
u - z ' ' -0.012 


_1 lfi+ 011 
1 - lu -0.11 


0.718 


0.0222 


3.8 


SNLS3+BAO+WMAP7+H0 


u.z i y_ .oo9 


1 1 9+O.O8 
x - "-0.07 


0.713 


0.0223 


425.1 


SNLS3+BAO+WMAP7+H0+Growth 


0-2801°;^ 


_n qq+ 005 


0.697 


0.0225 


440.2 



Notes. SNLS3 with systematics included, a and /3 fixed to their marginalized value, h and Q®h 2 have been minimized so no error 
bars are provided. 



5. Discussion 

In this section we compare our results to other recent pub- 
lications on the same subject. 

lApplebv fc Linden (|2012"D conclude that the uncoupled 
Galileon model is ruled out by current data since their best 
fit yields A\ 2 = 31 with respect to the best fit ACDM 
model. In addition, they obtain a long narrow region of 
degenerate scenarios with nearly the same likelihood. In 
our case, the best fit has A% 2 = 10.2, we obtain enclosed 
contours in all projections and a clear minimum. 

Despite the fact that we are using the same ex p ansio n 
and perturbation equations as lApplebv fc Linden (|2012fl . 
there are differences between the two works. We use a 
parametrization of the model which makes our study in- 
dependent of initial conditions for x, while they set X4 — 
x(zi = 10 6 ) by imposing a p-n{zi) that varies in their pa- 
rameter scan. This requires one to solve a S'^-order polyno- 
mial equation in Xi - and hence to be forced to choose one 
of the five solutions - or to assume one of the four terms 
CiH 2+2 ( l ~ 2 ) is dominant in the (00) Einstein equation. 
In any case, this leads to a parameter space which is dif- 
ferent than the one we explored. Another difference arises 
from the theoretical constraints which are used to restrict 
the parameter space to viable scenarios only. Our set of 
theoretical constraints is larger since we also use tensorial 
constraints, which proved to be very powerful. This also 
leads to a different explored parameter space. 

In iDe Felice fc Tsuiikawal (|2010D . the rescaling of the 
Galileon parameters is done using a de Sitter solution in- 
stead of using xo as in this paper. This leads to relations 
fixing their "62" and "C3" coefficients as a function of their 
"cV and "C5" coefficients (denoted a and /3 in their study), 
but requires two initial conditions to compute the cosmo- 
logical evolution. Those are also fitted using experimental 
data. W ith this param e trizat ion and without growth con- 
straints, iNesseris et al.l (|2010f) find best fit values for their 
"C4" and "C5" of the same sign and same order of magni- 
tude as in our work, despite our d i fferent paramctrizations. 
A second paper by lOkada et al.l ((2012) includes redshift 
space distortion measurements and rules out the Galileon 
model at the lOcr level. 

The first difference with respect to our work is in the 
treatment of the initial conditions, and the use of an ex- 



tra theoretical constraint to avoid numerical instabilities 
during the transition from the matter era to the de Sitter 
epoch. This reduces the parameter space with respect to 
that explored in our work. As already stated above, we 
think that a better modelling of G^g including non linear 
effects should be conducted instead of discarding scenarios 

with such in stabilities. 

Second, lOkada et al.l (|2012j) use fas measurements 
not corrected for the Alcock-Paczynski effect. Moreover, 
to make th e ir fa s predictions in the Galileon model, 
lOkada et al.l (|2012t ) sets the normalisation of ag today to 
the WMAP7 <j$(z = 0) measurement, which was obtained 
in a cosmological fit to the ACDM model. This normalisa- 
tion leads to the following <j%{z) evolution: 

8 i z ) = a S ( 'iXO) ' ' 

This assumes that the Galileon theory predicts a matter 
power spectrum s imilar to that of A CDM at z = 0, which 
is not guaranteed (jBarreia et al J . 120121) . On the contrary, we 
use the WMAP7 a% measurement to set the normalisation 
at decoupling z ss z* (see Equation |55|) . Thus we take into 
account the different growth histories between the ACDM 
and the Galileon models (Equation [60] which is different 
from Equation I68[) . We can compare our best fit scenar- 
ios for these two models to the fas and F measurements. 
Figures [12] and [13] show the result of this comparison. The 
agreement with data is good in both models. In particu- 
larly, the Galileon model does n ot exhibit a d i screp ancy as 
strong as was found in Fig. 3 of lOkada et al.l ()2012l ). 



6. Conclusion 

In this paper we have confronted the uncoupled Galileon 
model with the most recent cosmological data. We intro- 
duced a renormalization of the Galileon parameters by the 
derivative of the Galileon field normalized to the Planck 
mass, in order to break some degeneracies inherent to the 
model. Theoretical conditions were added to restrict the 
analysis to viable scenarios only. This allowed us to break 
the parameter degeneracies which otherwise would have 
prevented us from obtaining enclosed probability contours. 
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Fig. 12. fas(z) measurements from different surveys (6dF- 
GRS, 2fFGRS, SDSS LRG, BOSS and WiggleZ) compared 
with predictions for the ACDM model (with parameters of 
Table [5]- dashed purple line) Galileon scenarios. The solid 
blue line stands for the best fit Galileon scenario using all 
data, whereas the orange dashed line stands for the best fit 
Galileon scenario using growth data only. 
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Z 

Fig. 13. F(z) measurements from different surveys (BOSS 
and WiggleZ) compared with prediction for the ACDM 
model (with parameters of Table [5] - dashed purple line) 
and for Galileon scenarios. The solid blue line stands for 
the best fit Galileon scenario using all data, whereas the 
orange dashed line stands for the best fit Galileon scenario 
using growth data only. 



In particular, the conditions on the tensorial propagation 
mode of the perturbed metric proved to be very helpful. 

We used a grid search technique to explore the Galileon 
parameter space. Our data set encompassed the SNLS3 
SN la sample, WMAP7 {Z a ,i?,2:*} constraints, BAO 
measurements and growth data with the Alcock-Paczynski 
effect taken into account. We found {£^,02,03,04} = 

In 97-1+O.OI3 a oco+0.518 1 rn 7 +0.203 n771+°- 098 \ 
\ u - z ' i -0.008' *- ooz -1.220> 1 ' OB ' -0.726' u - ' ' i -0.061 J • 

The final x 2 is slightly above that of the ACDM model 
due to a poorer fit to growth data. 

The best fit Galileon scenario mimics a ACDM model 
with the three periods of radiation, matter and dark en- 



ergy domination, with an evolving dark energy equation of 
state parameter w(z), and an effective gravitational cou- 
pling G^g(z). Predictions for the latter are possible only 
in the linear regime, which may have an impact on our 
results derived from growth data as the latter were com- 
puted using a non-linear theory. A more precise theoretical 
and phenomenological study should be conducted to fairly 
compare the Galileon model with these data. 

Our best fit is more favourable to the Galileon model 
than other recent results. The main difference between our 
treatment and those works lies in the treatment of initial 
conditions. We have also tried to make as few assumptions 
and approximations as possible when computing observ- 
able quantities. Finally, when using growth data, we paid 
attention to choose measurements which were derived in a 
model-independent way. In the future, a study considering 
pr ecise predictions of t he full power spectra as suggested 
bv lBarreia et all (|2012f) would provide more stringent tests 
of the validity of the Galileon model. 
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Appendix A: Instability of probability contours 

Instead of absorbing the initial condition xq in the Ci — > 
Ci redefinition, we can be tempted to fix it using the 
(00) Einstein equation at z=0 for each scenario: 



1 - fi° - - 



-2c 3 ^-^fc 4 4 + 7c5X 5 = (A.l) 



1 a „ 15 
gC 2 x + 2c 3 x - — 

To find xq, a fifth-order polynomial equation is to be 
solved, which can lead to 5 complex solutions. A reasonable 
choice is to keep only the scenarios that give a unique real 
solution. 

The system of differential equations [3] and |H] adopts an 
unusual behaviour. Referring to Fig. IA.11 the shape of the 
probability contours remains unchanged whatever the lim- 
its of the scanned parameter space. In other words, the 
likelihood surface is invariant when the limits of the ex- 
plored parameter space are proportionally changed. The 
model seems to exhibit a scale invariance allowing data to 
be fitted whatever the boundaries of the explored parameter 
range. Moreover, we cannot obtain contours well enclosed 
in any explored parameter space: the likelihood surface has 
an infinite valley of minimum \ 2 instead of a unique mini- 
mum. 

Equation IA.1I shows that small c^s produce a high xq, 
and high CjS a low x$. Nevertheless, the theoretical con- 
straints of Sect. 12.51 cannot favour or disfavour high CjS or 
Xq as they contain also this correspondence between the CjS 
and Xq. So, for different sets of CjS, identical cosmological 
scenarios are computed whatever the scale of the c^s: the 
important point is that these equivalent scenarios have the 
same Q„(z), whether this is due to high or small qs and as 
a consequence have the same H(z) evolution and then the 
same x 2 - 

Thus, a scale choice has to be done to fix the likelihood 
surface, but this choice has not to be arbitrary. A solution 
is provided in 12. 21 by absorbing Xq into new parameters CjS. 
This new parametrization absorb a degree of freedom and 
allow us to use (00) Einstein equation to fix C5. This may 
b e the origin of the degene racy in \ 2 reported in Sect. Ill 
of lAnnlehv fc Linderl (|2012h . 
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Fig. A.l. Experimental constraints on the Galileon model 
from SNLS3 data for different ranges in as using the 
method developed in Appendix A, with a and fix ed to 
their ACDM best fit values from lSullivan et all ()201ll ). The 
four-dimensional likelihood £($7^,02,03,04) (05 fixed to 
here) is marginalized over C3, C4, .A/J^, M 2 B to visualize the 
O^i, C2 contour plots. The filled dark, medium and light 
blue contours enclose 68.3, 95.4 and 99.7% of the probabil- 
ity, respectively. 



Appendix B: Approximation for /„ computation 

The computation of l a (see equation 1421) requires the evo- 
lution of the cosmological model from today to a — (see 
equation HQ]). In the Galileon context, the non-linear evo- 
lution equations require more and more precision and finer 
steps when approaching the limit a — > 0. In addition, it is 
physically questionable to extrapolate the Galileon model 
up to the very first instants of the Universe. 

So our iterative computation is stopped at a certain 
a-min close to 0, without affecting significantly the final 
value of I a - Let a m i„_i be the step before a m i n where the 
cosmological equations are computed, and f(a) the inte- 
grand function of r s (z»). Although the integral is stopped 
at a = a rn i„, we can compensate this approximation by a 
linear interpolation of the integral : 
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Fig. A. 2. Left panel: evolution of l a with a m i n without the linear interpolation as described in the text, for a subset of 
Galileon scenarios. Note that most scenarios approach the WMAP7 measurement l a ss 300. Right panel: correction to l a 
for different values of a m i n and for the same subset of scenarios as in the left panel. The dashed line is the value of oi a , 
the WMAP7 measurement error on l a . 
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In the left panel of Fig. IA.21 we present the evolution 
of l a with a m in without the linear interpolation, for a sub- 
set of Galileon scenarios. The smooth evolution with a m i n 
allows us to consider the linear interpolation as a reason- 
able assumption. Moreover, for a m .;„ < 10~ 6 , the value 
of l a changes less than the WMAP7 measurement error 
<7; a = 0.76, as shown in the right panel of Fig. IA.21 Based 
on these results, we decide to use a m i n = 10 -7 , which pro- 
vides a correction on l a an order of magnitude below o\ a . 
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